A fast forward computational method for nuclear measurement using volumetric detection 
constraints 


Qiong Zhang", Lin LvLin' 
! University of Electronic Science and Technology of China, Chengdu 610000, China 
Corresponding author: zhanqio@uestc.edu.cn 


This work is supported by the China Natural Science Fund (U23B20151 and 52171253). 


Abstract 


Owing to the complex lithology of unconventional reservoirs, field interpreters usually need to provide 
a basis for interpretation using logging simulation models. Among the various detection tools that use 
nuclear sources, the detector response can reflect various types of information of the medium. The 
Monte Carlo method is one of the primary methods used to obtain nuclear detection responses in 
complex environments. However, this requires a computational process with extensive random 
sampling, consumes considerable resources, and does not provide real-time response results. Therefore, 
a novel fast forward computational method (FFCM) for nuclear measurement that uses volumetric 
detection constraints to rapidly calculate the detector response in various complex environments is 
proposed. First, the data library required for the FFCM is built by collecting the detection volume, 
detector counts, and flux sensitivity functions (FSF) through a Monte Carlo simulation. Then, based 
on perturbation theory and the Rytov approximation, a model for the detector response is derived using 
the flux sensitivity function method and a one-group diffusion model. The environmental perturbation 
is constrained to optimize the model according to the tool structure and the impact of the formation 
and borehole within the effective detection volume. Finally, the method is applied to a neutron porosity 
tool for verification. In various complex simulation environments, the maximum relative error between 
the calculated porosity results of Monte Carlo and FFCM was 6.8096, with a root mean square error of 
0.62 p.u. In field well applications, the formation porosity model obtained using FFCM was in good 
agreement with the model obtained by interpreters, which demonstrates the validity and accuracy of 
the proposed method. 
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1. Introduction 


As the exploration of oil and gas reservoirs has gradually deepened, the use of unconventional 
reservoirs has been increasing in proportion and has gradually become one of the main resources for 
petroleum production. Nuclear measurements, typically obtained using neutron porosity and gamma 
density tools, play a critical role in reservoir evaluation. However, in unconventional reservoirs 
characterized by various lithologies, complex pore structures, diverse fluid types, heterogeneity, and 
thin layers [1,2], nuclear responses are heavily impacted, causing difficulty in obtaining accurate 
formation information from a single type of measurement [3]. To obtain reliable formation parameters, 
multiple logging tools based on different physical principles and comprehensive analysis are usually 
needed. Hu developed a multiphysics model to extract key features of logs, including natural gamma, 


density, neutrons, resistivity, and acoustics, to predict rock mechanical parameters [4]. Chen 
introduced a symbolic regression model that used different sources of data to evaluate reservoir 
heterogeneity [5]. The main drawback of these methods is that they require the use of multiple tools, 
which means that the measurement process can be expensive and heavily relies on the knowledge of 
tool operators. Therefore, an idea was proposed to determine the formation parameters by updating 
them in the simulation model until their response matches that of the measurement. Tang obtained the 
net inelastic and capture spectra based on Geant4 simulations to guide the separation of inelastic and 
capture gamma rays for measurements [6]. Mehana integrated production analysis with simulations to 
predict reservoir oil recovery [7]. Although these attempts are meaningful, obtaining reliable results 
from simulations in different logging environments can be extremely time-consuming, making it 
almost impossible to obtain nuclear tool responses in real time using conventional methods such as 
Monte Carlo [8]. 


In recent years, a fast forward computational method (FFCM) has been proposed. To realize real-time 
updates of the formation model, Shao developed a gamma-log-based FFCM and applied it in different 
downhole environments [9]. Qin proposed an FFCM by constructing a series of nonlinear functions to 
describe the logging responses for different boundary distances of formations [10,11]. Similarly, 
FFCM can also be used for the dynamic monitoring of directional wells [12]. Luycx and a research 
team from UT Austin introduced the concept of the flux scattering function into FFCM, which utilizes 
the spatial scattering flux to improve the accuracy of the simulated tool responses [13-15]. Liu 
proposed an improved flux sensitivity function (FSF) method by extending the Taylor expansion to 
the second order to minimize the perturbation of environmental changes, thus significantly improving 
accuracy [16]. However, the existing FFCMs do not consider the detection volume of detectors and 
only correct the FSFs of different detectors in an iterative manner. Therefore, by considering a priori 
information, that is, the detection volume of different detectors, a novel FFCM using volumetric 
detection constraints 1s proposed. The detection volume can be considered the maximum region that a 
detector can detect and can typically be defined by the depth of investigation pertaining to each detector. 
Because the detection volume contains key information, such as the direction of tool measurement and 
the angle ofthe wells, this enables the proposed method to accommodate more environmental variables 
simultaneously, and could be feasible in more complex environments. 


Section 2 presents relevant nuclear physics theories based on which an FFCM is developed. Section 3 
describes how the method is implemented in an existing neutron porosity tool and tested in both 
simulated and experimental environments. The results are discussed in detail. Section 4 presents the 
conclusions. 


2. Methodology Development 


In this section, an FFCM using volumetric detection constraints is designed; its implementation 
workflow is illustrated in Fig. 1. A preset data library, including the detection volume and FSF, is first 
built using Geant4 [17-21], a software that simulates the passage of particles through matter. A 
weighted migration length is obtained by inputting the environmental parameters of the target model, 
including lithology, borehole parameters, rock layer, and angle of deviated well. Then, the most 
suitable base case, including the corresponding correction function and flux sensitivity function FSF; 
are selected by the migration length using the distance method [22]. The weighted migration length is 
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input into the correction function to obtain the final corrected migration length. Finally, the detector 
counts are calculated by substituting the final migration length and FSF; into the detector response 
model, and the formation porosity is calculated based on the near-to-far detector count ratio. 
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Fig. 1 Workflow of FFCM constrained by detection volume implementation 
2.1. Detector response model development 


To obtain the model form of the detector response and ensure that the accuracy of the FFCM meets the 
requirements, a numerical model of the detector response was obtained based on the perturbation 
theory [23]. Perturbation theory is a mathematical method for determining the approximate solution of 
a problem under a small perturbation based on the exact solution of an existing relevant problem. Its 
underlying formula can be written as 


Ao 
AN 5, 


Ao 
Np = N AN(—) = Nz | 1 + ——— 1 
P Bt =- B + N; (1) 


where Np, Ng are the detector counts under perturbation P and base case B, respectively. AN is 
the change in counts caused by the environmental perturbation Ao. According to perturbation theory, 


as AN 0, Ao > 0 and N' = lim To- 


According to the Rytov approximation [24], Eq. (1) can be expressed as 
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(2) 


Np = Np'e 


Eq. (2) is expanded into a Taylor series at N = Ng (corresponding to the environmental parameter 
Og) as Shown in Eg. (3), with the expansion term o. The first-order Taylor expansion is given by Eq. 


(4). 
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N' in Eq. (4) is obtained based on nuclear physics. 


The detector count N can be expressed as the sum of the particles in space from the source weighted 
by their contribution to the detector count, as shown in Eg. (5) [25-27]. 


N = [ar | az (r, E,A|r3)D* (r, E, O|rg) = | ar | az Fsr (5) 


where © and ©* are the forward flux and adjoint flux, respectively. rs and rg are the source and 
adjoint source locations, respectively. The flux sensitivity function is defined as FSF = 
d(r, E, Alr 6 * (r, E, Q|rg). 


Based on perturbation theory, the FSF of the perturbation model is approximately equal to that of the 
base case, as shown in Eq. (6). 


N' = FSF x FSF, (6) 


According to Eqs. (4) and (6), the difference between the base and perturbed counts at a certain spatial 
position due to the single neutron energy is given by 


AN» singlePosition = N'— = FSF, — (7) 


Because the energy of the americium-beryllium neutron source is a fixed energy spectrum, it can be 
considered as a linear superposition of neutron sources with multiple single energies. The detector has 
a certain volume; therefore, it must be integrated over space, as shown in Eq. (8). 


Ao 
ANp = | ar | az AN» singlePosition = | ar | dE Po (8) 
According to Eqs. (4) and (8), the detector count in the case of perturbation is given by Eq. (9). 


Ao Ao 
Np = Np + | dr [ dE ror = Np + | ar | dE ost. E, Ars) 080, E, Arg) — (9) 
B B 


Ac . : . ; ' a "-— ; 
where = 1s derived based on the volumetric detection constraints method, which is described in detail 
B 


in Section 2.3; the base case's forward flux ©p and adjoint flux ©} are solved based on the time- 
independent Boltzmann eguation, as shown in Eg. (10) [28,29], which describes the entire process of 
neutrons attenuating in the formation and finally reaching the detector under neutron porosity logging 
using a radioactive neutron source. 


[Q V + u(r, E)]b (r, E, Q) = [ar f aw us(r,E' > E,Q > Q)O(r, E, Q0) + § (10) 


© is the neutron flux at position r in the angular direction Q and with energy E; p is the total 
reaction cross section, which can be expressed as the inverse of the migration length Ly; 
us (r, E' > E, Q' > O) is the scattering cross section from energy E" to E and angular direction from 
Q' to Q; and V is a Laplacian operator. 


Considering only the monoenergetic neutron group, P(r, E, Q) can be written as ®(r, Q). Based on 
the simple diffusion approximation, i.e., no dependence ofthe angular flux on the direction 0, Eq. (10) 
can be simplified to a neutron diffusion equation, as shown in Eq. (11) [30]. 


[DV? — u]ó(r) = -S (11) 
1 


D = — 12 
3h (12) 


where D is the diffusion coefficient, given by Eg. (12). 


Based on the one-group neutron diffusion model [31] and a neutron source located at rs, Eg. (11) can 
be expressed as 


[DV? — o] (r|r;) = -S (13) 


where o is the perturbation parameter with G = 1/L,, [32]. ®(r|rs) is the flux field existing at r 
generated by a positive point source of S intensity located at rz, i.e., the spatially diffuse scalar flux 
at observation position r. 


The neutron flux equation can be solved using Eq. (13) as expressed in Eq. (14). 
ó(r|r;) = — Sa -air-rsi (14) 
Ss? AnD |r — rs| 


Combining Eqs. (5) and (14), when the position and energy are fixed, the detector count N is related 
to the perturbation G = 1/L,, of the environment. 


The numerical model of the detector count can be expressed as 
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where S? and are the source and adjoint source of base case, respectively. LP, is the migration 
length of the standard model. LB,/AL,, is derived based on the detection volumetric constraints 
method, which is described in detail in next section. 


22. Detection volumetric constraints method 


The detection volume contains important information such as the tool position, tool direction, and 
angle of the wells. For example, in the case of a high-angle well, the detection volume of the tool is 
inclined and the proportion of each layer in the detection volume differs from that of vertical wells. 
Therefore, the detection volume can be introduced to constrain the effect of perturbation to obtain 


IB, /ALm in Eg. (15) and optimize the detector response model to achieve a high-accuracy calculation. 


An open-hole well consists mainly of two parts: a borehole filled with water and a formation. The 
hydrogen content of the formation rock mainly originates from the water in the pores; therefore, the 
hydrogen content of the borehole fluid is much larger than that of the formation rock, and the borehole 
fluid has a greater ability to slow down neutrons. Therefore, the particle transport path of neutrons in 
the borehole is shorter than that in the formation. According to Eq. (14), it can also be proved that the 
particle transport trace in the borehole is shorter when the neutron flux decays to the same value. This 
implies that the detection volume can be approximated as an ellipsoid, which is further demonstrated 
using the actual tool in Section 3.1. 


Assuming that the well slope angles in the X, Y, and Z directions are 0%, 0, and 6,, respectively, the 
detection volume equation is given by Eq. (16). 


V = (x —C 73 P 6-6) A = C)? (16) 
Detection — NC 4 CUERO c? 
(x', y, z') 
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Where Vpetection is the detection volume; (C, Cy, Cz) are the center coordinates of the detection 


volume; (x',y',z') are the coordinates of the detection volume equation after rotation; a, b, and c 
denote half of the length of the ellipsoid along the three coordinate axes. 


When the borehole diameter and lithology are fixed, the variation in counts is caused by formation 
porosity. Therefore, the corrected migration length value Lin ing: and the correction value ALm can 
be obtained using Eqs. (18) and (19), respectively. Owing to the difference in the mechanical structure, 
the correction values are different for various detectors of different tools. 


Np — Ng 


Ln fi = 
m final i dr f dE FSF, + Op (18) 


EV formation 
Alm = Lm final = Ln (19) 


where I2, is the initial migration length calculated in an infinite medium with a point source. 


Based on Eg. (19), the migration length correction formula f for various detectors in each base case 
can be fitted using Eq. (20), and the fitted data are the eigenvalues and adjacent data of the standard 
model. 


Lm. final = Ln + ALm = f (Lm) (20) 


When the vertical resolution of the logging tool is not sufficiently high, or when thin layers are present, 
the tool detects signals from multiple geological layers simultaneously, as shown in Fig. 2, which is 


commonly referred to as a pile-up. This may result in logging data that do not accurately represent the 
characteristics of individual layers but rather a combination of their characteristics. Therefore, to 
accurately calculate the detector responses, it is necessary to consider the influence weight of the 
different geological layers detected by the detector in the response, which is achieved by weighting 
the migration length. Within the effective detection volume, let the detected volume of the i^^ rock 
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layer be Veo. ; (i = 1,2,--,n), the corresponding migration length of the it rock layer be 


La — (i = 1,2, ---, n), and the detected volume of the borehole liquid be Vzj,;4. A two-dimensional 


schematic of the detection volume for multiple rock layers is shown in Fig. 2. 


Fig. 2 Two-dimensional detection volume schematic 


The detection volumes for different rock layers and borehole fluids can be calculated using Eqs. (21) 
and (22), respectively. 
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BD 
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where BD denotes the borehole diameter and N denotes the intersection of the detection volume 
with the formation or borehole. 


When the borehole diameter is not equal to the standard borehole diameter, the borehole fluid volume 
around the tool changes, affecting the neutron transport and causing a variation in the detector count. 
The standard model is based on a fixed borehole diameter and does not consider the influence of the 
borehole fluid. When the borehole diameter changes, the weight of the borehole fluid influence is no 
longer equal to zero. A migration length weighting method is given by Eq. (23). 


Lin = WF luid i (Linor m Uo rua) + Uso; (23) 
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where Ll, is the weighted migration length of the impact from the borehole diameter constrained by 
the detection volume. wftuia is the weight reflecting the variation of borehole diameter and can be 
calculated using Eq. (24). The function g can be obtained by fitting the relative variation of the 


borehole fluid volume and the formation migration length. (L9,, — L? ) represents the effect of 


Mfor Th fluid 
borehole fluid and formation replacing each other. 


When the detector detects multiple layers of rock, the migration length must be weighted based on the 
detection volume of the different rock layers and L1,. The weighting equation is given by Eq. (25). 
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where 12, is the weighted migration length of multiformation impact. Lhi is the weighted migration 
length of borehole diameter impact of the i^" formation. Wfor_i 1$ the weight of the it formation. 


Combining Eqs. (20), (23), and (25), the final migration length can be obtained as 
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To select the base case that best matches the environmental perturbation, the distance method shown 
in Eq. (28) is used. 


Basecase = arg min (L5, — Lm final) (28) 


According to Eqs. (9) and (27), the detector response using volumetric detection constraints is given 
by 
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3. Verification 


3.1. Implementation in a neutron porosity tool 


Neutron porosity logging is an important technique for evaluating reservoirs. This method is applicable 
to any wireline porosity tool, and is discussed in this study based on a wireline neutron porosity tool. 
The neutron porosity tool was offset to obtain more information about the formation during logging. 
As shown in Fig. 3, the tool with a **'Am-Be source is placed in a fluid-filled borehole and attached 
to the borehole wall [33]. Two thermal neutron detectors filled with "He are placed 380.2 mm and 618 
mm from the source, and a shield made of a tungsten-nickel alloy is placed between the source and 
detectors to shield the axial gamma rays. The detailed parameters of the neutron porosity tool are listed 
in Table 1. During logging, neutrons are emitted into the formation and reach the detector for collection 
through processes such as elastic scattering, inelastic scattering, and capture. To simulate deep 
penetration problems, the Geant4 physics list called "shielding" is used, which is specifically designed 
for this purpose. 


Table 1 Neutron porosity tool parameters 


Component Materials Density ! m Distance to source 
(diameter*height) 
Near detector ?He 1.50 g/cm? 26.0mm*104.3nm 380.2 mm 
Far detector ?He 1.50 g/cm? 44.8mm*237.6mm 617.6 mm 
Detector sleeve Steel 4.50 g/cm? 56.5 mm*651.0 mm 634.5 mm 
Detector housing Steel 7.93 g/cm? 73.0 mm*944.0 mm 746.2 mm 
Shield Tungsten-nickel alloy 18.40 g/cm? 34.1 mm*24.4 mm 157.1 mm 
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Fig. 3 Neutron porosity logging model 


To ensure the accuracy of the model in a short computation time, the neutron logging model is 
described by a two-dimensional Cartesian coordinate system [34], that is, integration over X, where 
the formation is vertical to the Z axis and the neutron porosity tool is symmetrical about the Z-Y plane. 
The space around the tool is divided into grids with axial and radial boundaries to collect neutron flux. 
To build an FFCM covering multiple environments, a data library containing various parameters of 
borehole diameter, borehole fluid salinity, and porosity was constructed, as listed in Table 2. The 
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corresponding migration length values were calculated for different formation porosities. In this case, 
the standard borehole diameter for the base case was 311.2 mm and other environmental parameters 
ofthe base case are listed in Table 2. The total computational time for the database was 210 h. It should 
be noted that as the number of environmental factors increases, larger databases are required to ensure 
computational accuracy. 


Table 2 Database parameters for the FFCM 
Lithology Borehole diameter (mm) Salinity (kppm) Porosity (p.u.) 
limestone 215.9, 431.8, 311.2 0, 100, 200, 300 0.9, 5, 10, 20, 30, 40, 50, 60 


In this study, Geant4 was used to construct a neutron porosity logging model and obtain a data library. 
The conversion equation for the detector ratio and porosity is given by Eq. (30) [35,36] and can be 
obtained from the detector ratio data in the data library. 


bp = 1.16: R? + 2.05- R — 8.11 30) 


where $ is the measured porosity and R is the near-to-far detector ratio. 


The established data library contains the response of the model in different environments, including 
the forward flux, adjoint flux, and FSFs. The two-dimensional distributions of the forward and adjoint 
fluxes were obtained using Eq. (14) [37]. To better show the distribution, the normalized flux 
distribution and FSF are shown in Figs. 4-8. As shown in Fig. 4, the forward neutron flux gradually 
decreases in all directions centered at the location of the source, and is more widely distributed in the 
formation direction because the slowing ability of the formation to neutrons is weaker than that of the 
borehole fluid. The adjoint flux is acquired by setting the adjoint source on the detector; thus, the 
distribution of the flux centered on the detector and decreasing in all directions is presented in Figs. 5 
and 6. Because the increased porosity of the formation increases the hydrogen content of the formation 
and thus enhances the scattering effect on neutrons, regardless of the forward flux or adjoint flux, the 
normalized neutron flux distribution is more concentrated at the source center or adjoint source center 
with increasing porosity of the formation. As shown in Fig. 7 and 8, greater flux sensitivity values are 
available in the source and detector regions because these regions have the highest neutron flux 
recorded in the detector, which reflects the contribution of neutrons to the detector in each region. The 
flux sensitivity functions reflect the common characteristics of the forward and adjoint fluxes; that 1s, 
the sensitivity values are more concentrated near the source and detector as the porosity of the 
formation increases, and the sensitivity values are not concentrated in the part of the tool between the 
source and detector because of the shielding effect inside the tool [38]. 
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Fig. 4 Normalized forward flux distribution at 311.2 mm borehole diameter with different formation porosity 
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Fig. 5 Normalized adjoint flux distribution of Near Detector at 311.2 mm borehole diameter with different 
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Fig. 6 Normalized adjoint flux distribution of Far Detector at 311.2 mm borehole diameter with different 
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Fig. 8 Normalized flux sensitivity function of Far Detector at 311.2 mm borehole diameter with different 
formation porosity 


After obtaining the FSF, it is necessary to obtain the three parameters, a, b,and c for the detection 
volume in Eq. (16). The transport tracks of the neutron particles that arrived at the detector were 
collected using Geant4. The maximum distance reached by 90% of the detected neutrons was defined 
as the maximum detection volume. The a, b, and c of various detectors in different environments 
are listed in Table 3. 


Table3 a, b and c of near and far detectors in different environments 


: Formation Near detector Far detector 
Borehole Environment . f 
n . migration 
diameter (material; tiui a b C a b C 
i en 
(mm) proportion) ex (mm) (mm) (mm) (mm) (mm) (mm) 
Limestone, 
99.9196: H20, 2592 36417 324.45 482.52 50515 41948 74101 
0.09% 
311.2 Limestone, 8094; 
iacob 1432 298.719 27034 44611 41336 362.691 680.68 
H2O, 20% 
Limestone, 60%; 
HDeSone, 075 1106 27189 26212 43441 38492 234383 6571 
H30, 40% 
431.8 Limestone, 80%; 14.32 2921 268.56 445.15 406.79 35612 67935 
215.9 H20, 20% 1432 30422 27221 448.67 420.58 36709 68447 


From Table 3, without considering the angle of the well and direction ofthe tool, the detection volume 
is mainly related to the location of the detector and formation migration length, which are slightly 
influenced by the borehole diameter because the tool is attached to the borehole wall. Therefore, the 
detection volume can be determined based on the migration length of the base case to constrain 
environmental disturbance. Thus, a preset data library is built, and the FFCM can be computed 
according to the pseudocode shown in APPENDIX A. 


3.2. Simulated data verification 


To ensure the validity of the FFCM, a test database containing variations in environmental parameters, 
such as borehole diameter, lithology, and porosity, was established. Setting the number of source 
particles to 2 x 10? so that the statistical error in detector response counts is less than 2% and 
converges [39], and each model requires more than 5 h of computing using 110 threads on Intel© 
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Xeon(R) Platinum 8280 CPU (© 2.60GHz x 224. The formation in the logging model is composed of 
alternating sandstone or limestone with different porosities, and the borehole is filled with fresh water. 
The wireline neutron tool moved continuously from bottom to top, sampling the formation at 100 mm 
intervals to obtain the detector response. The method was validated in two cases: a vertical well (Case 
1) and a high-angle well (Case 2). 


Case 1 was designed using a model with different lithologies, formation thicknesses, and porosities for 
a nonstandard borehole diameter of 350 mm. A schematic of the model is shown in Fig. 9, and its 
specific parameters are listed in Table 4. The calculation results of the FFCM in Case 1 are shown in 
Fig. 10, where Track 1 and Track 2 are the normalized counts of the near and far detectors, respectively, 
and the normalization formula is shown in Eq. (31), which represents the probability of a neutron 
emitted from the source reaching the detector. MC Near Det Counts and FFCM Near Det Counts 
represent the counts of near detectors from MC and FFCM, respectively, and MC Far Det Counts and 
FFCM Far Det Counts represent the counts of far detectors from MC and FFCM, respectively. Track 
3 shows the formation porosity logs of the MC and FFCM obtained using Eq. (30), which are called 
the MC Porosity and FFCM Porosity, respectively. Track 4 shows the absolute and relative errors of 
the FFCM and MC simulation results, where the red lines are the error lines of 1 p.u. and 8%, 
respectively. 


N 
Nnormalized = Ns G1) 


where Nyormalizea is the normalized count. Ng; and N are the number of source particles and the 
detector count, respectively. 


Fig. 9 Case 1 model schematic 


Table 4 The formation parameters in Case 1 
Layer Thickness(cm) Composition (material, mass fraction) L9, (cm) 
1 40 Limestone, 79%; H20, 21% 14.11 
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2 100 Sandstone, 84%; H20, 16% 15.41 
3 100 Limestone, 72%; H20, 28% 12.69 
4 40 Sandstone, 67%; H20, 33% 11.92 
5 80 Limestone, 79%; H20, 21% 14.11 
6 90 Limestone, 94%; H20, 6% 20.07 
7 40 Limestone, 81%; H20, 19% 14.54 


Track 2 Track 4 
MC Near Det Counts MC Far Det Counts Absolute Error 
[o slo 0.0001 


FFCM Poros Relative Error 
B © — O2). 


E— wo ——H 0 


Fig. 10 Calculation results of the FFCM in Case 1 


As shown in Fig. 10, 50 measurement points were obtained by sampling at a sampling interval of 100 
mm in the Case 1 formation model at a depth of 4.9 m. As shown in Tracks 1-3, the results of the 
FFCM calculations are in good agreement with the MC simulation results in terms of trends and values, 
both for detector counts and formation porosity. As shown in Track 4, the absolute error of FFCM is 
mostly below 1 p.u., and the relative error is mostly below 6%, with a maximum absolute error of 1.19 
p.u., the maximum relative error of 6.80%, and the root mean square error (RMSE) of 0.62 p.u. 
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Case 2 is a high-angle well model with the Y-axis angle 0, being 10°, while the other axis tilt angles 
0, and 8, are zero and borehole diameter is 311.2 mm. The working state of the wireline neutron log 
was simulated to study the results of FFCM with different formation thicknesses and porosities, as 
shown in Fig. 12. A schematic of the Case 2 model is presented in Fig. 11 and Table 5. 


Fig. 11 Case 2 model schematic 


Table 5 The formation parameters in Case 2 


Layer Thickness(cm) Composition (material, mass fraction) L9, (cm) 
1 20 Limestone, 88%; H20, 12% 16.83 
2 80 Limestone, 83%; H20, 17% 15.12 
3 140 Limestone, 93%; H20, 7% 19.47 
4 30 Limestone, 69%; H20, 31% 12.21 
5 80 Limestone, 82%; H20, 18% 14.86 
6 70 Limestone, 78%; H20, 22% 13.85 
7 30 Limestone, 86%; H2O, 14% 16.06 
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MC Near Det Counts MC Far Det Counts i Absolute Error 
0001 0.00025 |0 0.000 0 


0 pu 25 


FFCM Near Det Counts | FFCM Far Det Counts FFCM Porosi Relative Error 
(0.0001 0.00025|0 0 p % 


Fig. 12 Calculation results of the FFCM in Case 2 


As shown in Fig. 12, 46 measurement points were obtained by sampling at a sampling interval of 100 
mm in the Case 2 formation model, which had a depth of 4.5 m. In Tracks 1-3, the calculated logs of 
FFCM and MC coincide well and show a good matching performance. According to Track 4, the 
absolute error of FFCM is mostly below 1 p.u. and the relative error is mostly below 6%, with a 
maximum absolute error of 1.09 p.u., a maximum relative error of 6.54%, and an RMSE of 0.56 p.u. 
Based on the validation of two simulation cases, the good agreement between FFCM and simulations 
verifies the validity and reliability of the novel algorithm. 


Once the database has been constructed, FFCM takes less than 1 s to execute, which provides a 
significant advantage in complicated scenarios that typically require heavy amounts of Monte Carlo 
modeling. For instance, in Case 2, the conventional Monte Carlo method requires the computation of 
a minimum of 46 models with the correct parameters, which would take 230 h with 110 threads on 
Intel® Xeon(R) Platinum 8280 CPU @ 2.60GHz x 224. Considering the additional calculations of 
other model parameters during the curve-matching process or longer, well-segmented processing 
requirements, the Monte Carlo method requires more time. However, using FFCM, each model in 
Cases 1 and 2 can be computed in less than 1 s regardless of the computational hardware parameters, 
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enabling real-time data processing. 
3.3. Experimental part: field application 


To evaluate the performance of FFCM, it was applied to a field well. The experimental well is located 
in the southern part of Liaodong Bay in the Bohai oil field and is drilled with a 241.3 mm (9.5 inch) 
bit. During logging, the tool slides against the borehole wall from the bottom to the top and is sampled 
at 76 mm intervals. The sieve residue log shows that the logging section is dominated by sandy 
mudstone formations, and borehole fluid sampling shows that borehole fluid has a density of 1.18 
g/cm? with a salinity of 94 kppm. The compositions of the formations with different porosities and 
borehole fluids were obtained to calculate the initial migration length L9,. In addition, the 
experimental well requires real-time measurements of the borehole diameter and angle of the well, 
which provides more accurate inputs to the method. The pseudocode of the method is shown in 
APPENDIX A, and its input is a formation model that includes the lithology, borehole parameters, and 
formation properties. By iteratively updating the input formation model to ensure that the calculated 
porosity log of the FFCM matches the measured porosity log of the tool as much as possible, the final 
porosity model with multiple layers can be used to assist in the rapid completion of formation analysis 
to support the formation of a layered interpretation. By utilizing the FFCM, as described above, the 
limitations of the baseline calculation method can be overcome, and the pile-up problem can be 
resolved, resulting in accurate properties of different geological layers. However, although the 
proposed method may introduce some errors compared to a full Monte Carlo simulation, it remains 
the optimal choice that balances time and accuracy considerations. 


The formation model obtained using the FFCM method in the X475-X515 well section was compared 
with the interpretation results, as shown in Fig. 13. Track 1 shows the borehole diameter log (CAL), 
natural gamma log (GR), and DEVOD log, reflecting the angle of the well. RD, RS, and RMSF in 
Track 2 denote the investigated deep double lateral resistivity log, shallow investigated double lateral 
resistivity log, and micro-spherically focused resistivity log, respectively. DT24 in Track 3 represents 
the acoustic slowness. Measured Porosity and Measured Density are the measured values from the 
neutron porosity and gamma density tools, respectively. Tracks 2 and 3 are important bases for the 
interpretation. The FFCM-calculated porosity (FFCM Porosity) log in Track 4 was obtained by 
iteratively and continuously updating the formation model, which matched well with the measured 
porosity. The FFCM Formation Model of Track 5 is the final formation model corresponding to Track 
4, and the Interpretation Model is the formation model obtained by synthesizing and analyzing a variety 
of measured logs; its detailed parameters are listed in Table 6. 


Table 6 The formation parameters from interpretation 


Laver Depth GR DT24 RD RS Density Porosity 
(m) (API) Cus/ft) (Q-m) (Q-m) (g/cm3) (p.u.) 
1 X473.0-X478.0 60.0 117.9 3.3 231 2.22 25.0 
2 X481.0-X482.0 65.7 109.5 3.8 33 22 22.6 
3 X487.9-X496.4 76.7 102.9 3.8 33 2.35 272 
4 X496.4-X499.4 58.5 100.7 16.0 14.3 2.27 24.9 
5 X499.4-X508.7 75.3 100.5 4.8 4.] 2.33 29.5 
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Fig. 13 Calculation results of the FFCM in field 


The GR log value of Track 1 in Fig. 13 was approximately 75 API, which is consistent with the sieve 
residue log, indicating that the lithology was mainly composed of sandy mudstone. The borehole size 
of the measurement section was maintained at 9—13 in because of the enlargement of the borehole, and 
the angle of the well in this section increased from 25°to 30°, as indicated by the DEVOD log of Track 
1. After multiple updates of the formation model, the FFCM-calculated porosity of Track 4 matched 
the measured porosity well. In addition, the formation porosity model of Track 5 was in good 
agreement with the model results provided by the interpreters, proving the validity and reliability of 
the FFCM constrained by the detection volume. 


4. Conclusions 


This paper proposes an FFCM that uses the detection volume to constrain environmental perturbation 
to optimize the detector response model, which is derived from perturbation theory and the Rytov 
approximation. The proposed FFCM was verified using a neutron porosity tool. Two cases were 
simulated based on various environmental disturbances, including the porosity, borehole, lithology, 
and slope angle of the well. The maximum relative error of porosity between FFCM and Monte Carlo 
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simulations was 6.80% and RMSE was 0.62 p.u. Regarding the field application, the formation model 
of FFCM matched with the interpretation results. The validity and reliability of the proposed FFCM 
were demonstrated. In addition, the proposed FFCM analyzed each model in less than 1 s, thus meeting 


the real-time application requirements for field use. 
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APPENDIX A 


Pseudocode of the FFCM using detection volumetric constraints for neutron porosity tool 


Algorithm 1: Fast forward computational model of neutron porosity based on detection volumetric 
constraints 


Input: lithology, distribution of rock layers and their porosity, borehole parameter, angle of deviated 
well. 
Output: counts of each detector, measurement porosity. 


1 


Calculate the initial migration length L9, based on the elemental composition of each 
formation and borehole fluid. 
Under the detection volumetric constraints, calculate the detected formation and borehole fluid 
volume according to Eq. (21) and Eq. (22). 
if (borehole diameter + standard borehole diameter) then 

Calculate the migration length LL, weighted by the influence of borehole diameter 
according to Eq. (23) and Eq. (24); 
else 


LL, =L 


Mor? 
end if 
if (detector detects multiple formation layers) then 
Calculate the migration length weighted by multiple rock layers according to equation Eq. 
(25); 
else 
pog 
end if 
Select base case through the distance method. 
Based on Eq. (27), the migration lengths Lm fina, of different detectors are calculate. 


Based on Eq. (29), the counts of different detectors and the near-to-far detector ratio are 
calculated. 
Based on Eq. (30), calculate the measurement porosity $y of the tool. 
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